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Abstract 

We demonstrate that simuhaneous reconstruction of scattering and absorption of a mesoscopic 
system using angularly-resolved measurements of scattered light intensity is possible. Image recon- 
struction is realized based on the algebraic inversion of a generalized Radon transform relating the 
scattering and absorption coefficients of the medium to the measured light intensity and derived 
using the single-scattering approximation to the radiative transport equation. 



I. INTRODUCTION 



There is considerable interest in the development of techniques for three-dimensional op- 
tical imaging of biological systems systems. In this context, of particular importance is 
the imaging of mesoscopic systems, where the photon transport mean free path is of the 
same order as the system size njl. For such systems none of the available optical ballis- 
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tic imaging modalities |2]-|6[ or the diffuse optical tomography [7| can be used. On the 
mesoscopic length scale, applications to biological systems include imagining of engineered 
tissues, semitransparent organisms, or superficial tissues. This article is the second in a 
series devoted to the problem of optical imaging in the mesoscopic scattering regime. In 
Ref . , we have proposed a novel imaging technique that uses angularly-selective measure- 
ments of scattered light intensity to reconstruct the attenuation coefficient of an inhomo- 
geneous medium, assuming that the single-light scattering is dominant. The advantages 
of this Single-Scattering Optical Tomography (SSOT) technique include the linearity, well- 
posedness, and two-dimensional character of the associated inverse problem, as well as the 
possibility to perform image reconstruction based on single-projection measurements. Here 
we generalize the SSOT technique to simultaneously reconstruct the internal scattering and 
absorption properties of the medium. 

We begin by presenting a brief review of the SSOT formalism introduced in Ref. 
We assume that the light transport in an inhomogeneous medium is described by the time- 
independent radiative transport equation (RTE) for the specific intensity /(r, s) of light at 
the position r and fiowing in the direction s, 

[s ■ V + /i„(r) + /i,(r)] J(r, s) = /i,(r) j A{s, s')/(r, s')d's' , reV . (1) 

Here /ia(i') and /is(r) are the absorption and scattering coefficients, and A{s, s') is the scat- 
tering kernels normalized such that J A{s,s')(Ps' = 1 for all s. The RTE is equivalent 
to the integral equation 
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where /;,(r, s) is the ballistic component of the specific intensity, and the ballistic Green's 
function Gb{r, s; r', s') is expressed as 
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the angularly- averaged ballistic Green's function. Here /if(r) = /ia(r)+/is(r) is the extinction 
(attenuation) coefficient. 

The light transport in a mesoscopic system is described by the ffist-order scattering 
approximation to the RTE. This corresponds to the assumption that light propagating in 
the inhomogeneous medium is just single scattered, and consists in replacing /(r, s) by 
/fe(r, s) in the right-hand side of Eq. ([2]). Consider that the medium is illuminated by a light 
beam of intensity Jq entering the slab at the point ri and in the direction Si, and that an 
angularly-selective detector registers the ray exiting the slab through the opposite surface 
at the point r2 and in the direction §2 (as shown in Fig. [T]). The intensity measured in a 
such experiment is denoted by /s(r2, §2; ri, §1). Within the single-scattering approximation, 
a relationship between the scattering and absorption coefficients of the medium and the 
measured light intensity is derived in the form 
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Here the integral Jggj^/ij(r(£))(i£ of the attenuation function is evaluated along the broken 
ray (BR) (shown in Fig. [1]), corresponding to single-scattered photons and uniquely defined 
by the source and detector positions and orientations, i is the linear coordinate on this ray, 
R21 is the ray turning point, and Jig is the average (background) value of the scattering 
coefficient. The data function 0(r2, §2; ri, Si) is defined as 

>2i sin 9i sin 6*2 / /s(r2, §2; ri, Si)(iv3s 
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where r2i = |r2 — ri|, the angles 9i and 62 are defined by cos6'i^2 = ^21 ■ Si,2, <^S2 is the polar 
angle of §2, and the scattering kernel A{s, s') is assumed position-independent and known. 

Eq. ([5]) is applied to optical imaging in the following manner. The selection of incidence 
and detection points and incidence and detection directions defines a slice in which image re- 
construction is performed. In Fig. [H this slice coincides with the FZ-plane of the laboratory 
frame. Assuming that the x-coordinate is fixed, the absorption and scattering coefficients 
can be regarded as two-dimensional functions of variables {y,z). On the other hand, the 
data are four-dimensional, depending, in general, on two spatial and two angular variables, 
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FIG. 1: (Color online) Geometrical illustration of the quantities used in Eq. ([6]) and elsewhere. 
"BR" denotes broken, single-scattered ray. 

corresponding to the source and detection y-coordinate, and source and detection direction, 
respectively. This enables the simultaneous reconstruction of scattering and absorption. 
By utilizing multiple incident beams and detecting light exiting the medium at different 
points, and by varying the incident and exit angles, it is possible to collect enough data to 
reconstruct the absorption and scattering coefficients in a given slice. Three-dimensional 
reconstruction is then performed slice-by-slice. 

In SSOT, simultaneous reconstruction of scattering and absorption can be in fact realized 
without scanning all parameter space. It is enough to keep the incidence direction (defined 
by the incidence angle Pi) fixed, to scan the incidence point yi, and for each such source 
realization to scan the detection point y2, for each detection position considering just one 
detection direction, corresponding to the detection angle (32 > Pi, if 2/2 > Ui + LtanjSi, or 
to the angle — /32, if 1/2 < yi + Ltan j3i, where (32 is fixed and L is the slab thickness. The 
incidence and detection angles are the angles between the 2;-axis of the laboratory frame 
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FIG. 2: (Color online) Schematic illustration of the proposed experiment geometry for simultaneous 
reconstruction of scattering and absorption. Reconstruction is performed in slices distributed along 
the X direction. The blue rectangles represent the areas in which reconstruction can be performed. 



and the unit vectors Si and §2, respectively. This principle is schematically illustrated in 
Fig. [2] for a rectangular sample illuminated by a normally incident beam. In the presence 
of scattering, both ballistic and scattered rays are present. To avoid the detection of the 
ballistic component of the transmitted light, the angularly-selective source and detectors are 
not aligned with each other. 

For a light beam entering the sample at position {yi, zi) and at the incident angle 
consider pairs of detections defined by the detection position and detection angle (1/2, ^2) 
and /32 > /3i, and {y'2, Z2) and —f32, respectively, where y2 > Vi + Ltan/3i, y'2 = y2 — 2{z2 — 
zq) tan P2, and Zq = {z2 tan P2 — Z1 tan /5i)/(tan /92 — tan Pi) — (1/2 — yi) /(tan /52 — tan This 
two measurements corresponds to rays single-scattered at the same position R21 = {yo, zq) 
within the sample, where yo = {z2 — zo)ta.ni32- Such rays are shown by the same color 
in Fig. O Writing Eq. ([5]) for the two situations and then subtracting the corresponding 
equations, one obtains the following equation: 

[ f,tW),z{i]di- [ fitW),z{i]di = (j){y2,P2;yi,Pi)-(p{y'2,-P2;yi,Pi] 

J BR(y2 Am ,/9i ) JBRiy'^-fS^wuM 

(7) 

Here the integrals of the total attenuation coefficient are along the two single-scattered 
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rays described above, and (j){y2, P2', Vi, Pi) and (j){y2, —P2',yi, Pi) are the data functions cor- 
responding to the two detections. Eq. ([7]) can be solved with respect to fit{y,z). Then 
Eq. ([5]) written for one of the detections can be solved for ln[/is(?/, z)/ (/is)]. From these two 
solutions, ^a{y, z) and Hsijl, z) can be determined. 

II. IMAGE RECONSTRUCTION 

In what follows, we illustrate simultaneous reconstruction of scattering and absorption in 
SSOT based on the algebraic inversion of Eq. ([5]). We note that a more sophisticated image 
reconstruction algorithm based on a inverse formula is also possible and will be presented 
elsewhere. 

The forward data is obtained by solving Eq. (I2j) numerically, along the lines presented in 
detail in Ref. j^, generalized here for the case of a homogeneously scattering medium. For 
isotropic scattering (yl(s,s) = l/47r), Eq. ([2]) takes the form 

/(r, s) = /,(r, s) + y G,(r, s; r', s')^^^u{v')d?r' cPs' . (8) 

Here u{y) = J J(r, s)d'^s is the density of electromagnetic energy satisfying the following 
integral equation 

n(r) = n,(r) + J g,{r, r')^u{r')d'r' , (9) 

where Uf,{r) = J If,{r,s)(Ps is the "ballistic density". The scattered component of the inten- 
sity. Is = I — lb, is computed by solving first Eq. ([9]) and then substituting the numerical 
solution u{r) into ([8]). We emphasize that this numerical approach is non-perturbative and 
includes all scattering orders, similarly to the experimental situation when all scattered light 
is detected. 

Equation is discretized on a rectangular grid and solved by methods of linear algebra. 
The energy density u{r) and the scattering and absorption coefficients are assumed constant 
within each cubic cell. The corresponding values m„ = M(r„), where r„ is the center of the 
n-th cubic cell, obey the algebraic system of equations 

(1 - RcqHsirn)) - 1— fl'6(rn, rm)fis{rm)Um = Mfe(r„) . (10) 

471 ^—f 
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Here h is the disctretization step, u;,(r„) = Jy Uf,{r)(Pr, and -Req = (3/47r)^/^/i is the 
radius of a sphere of equivalent volume to that of a cell, introduced to compute the diagonal 
matrix elements of the system (fTOjl 8|. The system of equations (11 01) is solved by direct 
matrix inversion, and then the specific intensity is calculated according with the discretized 
version of ([8]), 

/(r2,S2) = — gb{r2,rn)Hs{rn)un , (11) 

r2-r„=S2|r2-r„ 

where the summation is performed only over such cells that are intersected by the ray exiting 
from the detection point r2 in the direction §2. The data function is calculated in terms of 
the average of the specific intensity over the cell, /(r2, §2) = (l/h^) J dT2 /(r2, §2). To model 
noise in the measured data, /(r2,S2) was scaled and rounded off so that it was represented 
by 16-bit unsigned integers, similar to the measurement by digital ccd cameras. Then a 
statistically-independent positively-defined random variable was added to each measurement 
-^(r2,S2). The random variables were evenly distributed in the interval [0,nJav], where n is 
the noise level and Jav is the average measured intensity (a 16-bit integer). The date function 
is calculated using the discretized version of Eq. ([6]), 

'An /(r2,S2) 

Image reconstruction for the attenuation coefficient is obtained using Eqs. ([5]) and ([7j), 
which are discretized on the same grid as the one used for obtaining the forward solution, 
except that in this case planar slices with fixed x-coordinates are used. The discrete version 
of (HD is 

CunfJ'tn = 4>u , (13) 

n 

where the index z/ = (yi, A; ^2, /52) corresponds to a given realization of the source and 
detection pair, and C^n = — C^Jn, with the matrix element £in given by the length of 
the intersection of the detected ray i with the n-th cubic cell, (p^, = (f)l^^ — with (f)u^ 
the data function corresponding to the ray i. Eq. (ITB]) is solved for fim by regularized SVD 
pseudoinverse [9|], namely 

= iC*C)-'C\<f>) , (14) 

where 



(f>{y2,p2;yi,Pi) = - In 
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Here Q{x) is the step function, e is a small regularization parameter, and and (j„ are the 
singular functions and singular values, respectively, of the matrix C, obtained by solving the 
symmetric eigenproblem C*C\gn) = a1\gn)- Further, the scattering coefficient is determined 
from Eq. ([5]), discretized as 



Finally, the absorption coefficient is obtained as /lan = l^tn — l^sn- 
A. Numerical Results 

We considered a rectangular isotropically scattering sample of dimensions = 25h, 
Ly = 122h and = 40/i. The background scattering coefficient is chosen such that the 
optical depth fisLz is 1.6. This corresponds to the mesoscopic scattering regime in which the 
image reconstruction method of SSOT is applicable. The background absorption coefficient 
is set equal to fia = 0.1/2^. The target is a set of inclusions concentrated in the layers x = 6h, 
X = 13h and x = 20h. 

Image reconstruction is performed in slices x = XsUcc = const separated by the distance 
Ax = h. The reconstruction area inside each slice is 44h < y < 77h, Ah < z < 37h, with 
the field of view 34h x 34h. For each slice, the sources are normally incident {(3i = 0) on the 
surface ^ = at the positions x = XsHcc, y = Vs = nh, z = 0, with n integers. The detectors 
are placed on the opposite side of the sample at positions x = XsUcc, y = yd = nh, z = 
and measure the specific intensity exiting the surface 2; = at the angle of (32 = '?i"/4, for 
yd > ys, and — vr/4, for < yg, with respect to the 2;-axis. 

First, we considered the case of purely absorbing inhomogeneities, spatially modulated 
as shown in Fig. H] in the column marked "Model" . The scattering coefficient is constant 
throughout the sample and equal to the background value. The absorption coefficient for 
the inhomogeneities in the slice x = 6h is set to /la = 2fia and /i^ = 5/ia, for the outer 
and inner square, respectively. In slices x = 13h and x = 20h there are more absorbing 
inhomogeneities, of absorbing coefficient = 2/ia, 4/ia, 5/2^, going from the outmost to the 
innermost square. Thus, the contrast of fit (the ratio of fit in the target to the background 
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FIG. 3: (Color online) Image reconstruction for the total attenuation coefficient /it for a homoge- 
neously scattering sample with flsLz = 1-6 and fia = O.lfig, and for various noise levels n. The 
rows show the slices x = 6/i, 13h and 20/i, where the absorbing inhomogeneities are placed. The 
contrast in varies from 1.09 to 1.36. 

value) varies from 1.09 for the outmost squares in each shce to 1.36 for the innermost square. 
The results of image reconstruction for the total attenuation coefficient /xt and absorption 
coefficient fia for various noise levels n are presented in Figs. [3] and HI Only the slices 
containing inhomogeneities are shown. The other slices present no features, which means 
that no cross-talk exists between various slices, as it was also demonstrated previously j^. 
It can be seen that the spatial resolution of images depends on the noise level and can 
be as good as one discretization step, h. Note that image reconstruction is in very good 
quantitative agreement with the model (all panels in each figure are plotted using the same 
color scale) and stable in the presence of noise. Note also that the two-angle measurement 
scheme considered here enables better image reconstruction of the attenuation coefficient 
than the single-angle scheme considered previously sl. 

Consider now the case when scattering inhomogeneities are also present in the system and 
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FIG. 4: (Color online) Image reconstruction for the absorption coefficient /Za for a homogeneously 
scattering sample with /XsL^ = 1.6 and fia = O.l/Xs, and for various noise levels n. The rows show 
the slices x = 6h, 13h and 20h, where the absorbing inhomogeneities are placed. The contrast in 
Ha varies from 2 to 5. 

are spatially modulated as shown in Fig. [6] in the column marked "Model" , the absorption 
being modulated as described above. The scattering coefficients for the inhomogeneities in 
the slice x = 6h is set to fig = 1.33/is and fig = 1.66/2^, for the outer and inner square, 
respectively, and in this slice the the absorbing and scattering inhomogeneities overlap with 
each other. In the slice x = 13h, there are more scattering inhomogeneities as compared to 
the shce x = 6h, the scattering coefficient is fig = 1.33/2s, 1.66/is, 1.66fis, going from the out- 
ermost to the innermost inhomogeneity, and the absorbing and scattering inhomogeneities 
do not overlap. In the slice x = 20h, the absorbing and scattering inhomogeneities overlap, 
the scattering coefficient is modulated the same as in slice x = 13h except that its value 
for the innermost inhomogeneity is larger, fig = 2 fig. For this sample, the contrast of fit 
varied from 1.09 for the utmost squares in slice x = 13 to 2.27 for the innermost square in 
slice a; = 20. Imagine reconstruction in this case is presented in Figs. IMZl Very good image 
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FIG. 5: (Color online) Image reconstruction for the total attenuation coefficient /it for an inho- 
mogeneously scattering and absorbing sample and for various noise levels n. The rows show the 
slices X = 6h, 13h and 20h, where the inhomogeneities are placed. The background scattering and 
absorption coefficients are set such that flgLz = 1.6 and fia = O.l/Zs, the contrast in fig varies from 
1.33 to 2, the contrast in fia varies from 2 to 5, and the contrast in fj,t varies from 1.09 to 2.27. 

quality is obtained for both the total attenuation coefficient and scattering coefficient, image 
reconstruction for scattering being less influenced by the noise in the data function. On the 
other hand, image quality for absorption is notably lower. However, most of the relevant 
features are legible. 

Further, consider stronger scattering inhomogeneities and perform image reconstruction 
for a sample in which the scattering coefficient of the inhomogeneities is spatially modulated 
the same as in the previous case, but it is increased by a factor of 1.5, the absorbing 
inhomogeneities having the same characteristics. Thus, in this case, the contrast in the 
scattering coefficient varies from 2 for the outmost inhomogeneity in each slice to 3 for the 
innermost inhomogeneity in slice x = 20. The contrast in fit varies from 1.09 to 3.18. The 
results are presented in Figs. [HlfTOl Very good image reconstruction is obtained for both the 
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FIG. 6: (Color online) Image reconstruction for the scattering coefficient fig for an inhomogeneously 
scattering and absorbing sample and for various noise levels n. The rows show the slices x = 6h, 
13h and 20/i, where the inhomogeneities are placed. The background scattering and absorption 
coefficients are set such that fisLz = 1-6 and fia = O.l/is, the contrast in fig varies from 1.33 to 2, 
and the contrast in fj,a varies from 2 to 5. 

total attenuation and scattering coefficients, but imagine reconstruction for absorption is 
very poor. 

The reconstructed image quality is determined by two factors, the amount of scatter- 
ing in the system and the noise in the data. In particular, for stronger scattering, the 
single-scattering approximation we employ may be inaccurate, leading to poor image recon- 
struction. In order to separate the influence of these factors on the image quality, we perform 
image reconstruction based on a data function corresponding only to single-scattered light, 
obtained using the so-called inverse crime. This consists of generating data using the same 
model that the inverse solver is based on. Speciflcally, instead of solving RTE numerically 
and using the solution to calculate the data function according to the deflnition the data 
function is calculated from ([7]), derived within the single-scattering approximation of RTE, 
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FIG. 7: (Color online) image reconstruction for the absorption coefficient fXa for an inhomoge- 
neously scattering and absorbing sample and for various noise levels n. The rows show the slices 
X = 6h, 13h and 20/i, where the inhomogeneities are placed. The background scattering and ab- 
sorption coefficients are set such that fisLz = 1.6 and fla = O.lfig, the contrast in //^ varies from 
1.33 to 2, and the contrast in ^u^ varies from 2 to 5. 

by replacing the extinction and scattering coefficients by those of the model. In this case, the 
influence of the amount of scattering in the sample on the image quality is eliminated, the 
only influence coming from the noise in the data. Image reconstruction for the same sample 
that was analyzed in Figs. ISlfTU] is presented in Figs. [TTIfT^ By comparing these two sets of 
results, one can see that the imagine quality for the attenuation and scattering coefficients 
is almost the same in both cases. Therefore, it can be concluded that the single-scattering 
approximation works very well for the scattering strength considered and that, in this scat- 
tering regime, the most influence on the image quality comes from the noise in the data. 
Moreover, Figs. [TTlfT3l show that the reconstructed coefficients experience various levels of 
influence from the noise in the data. While the reconstructed attenuation and scattering 
coefficients are very little influenced by the noise, the absorption coefficient experiences a 
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FIG. 8: (Color online) Image reconstruction for the total attenuation coefficient /it for an inho- 
mogeneously scattering and absorbing sample and for various noise levels n. The rows show the 
slices X = 6h, 13h and 20h, where the inhomogeneities are placed. The background scattering and 
absorption coefficients are set such that flgLz = 1.6 and fia = O.l/Zs, the contrast in /i^ varies from 
2 to 3, the contrast in /i^ varies from 2 to 5, and the contrast in /i^ varies from 1.09 to 3.18. 

much stronger influence, imagine quality being very poor even for a noise level of 3% in the 
measured intensity. 

To explain the various levels of influence of the noise in the data on the image quality, 
we perform a rough error propagation analysis. Assume that the scattered intensity is 
determined with an error 6 Is- From Eq. f|T^ and the definition ([6]) of the data function, it 
follows that the noise in the data results in a noise Sfitn in the total attenuation coefficient 
given by 

Sf^tn = v/2 ^ f 5^ %^e(a^ - e)] . (17) 

Here, we have assumed that the relative error in determining the scattered intensity is the 
same for the symmetric measurements used to reconstruct /ij, and the factor results from 
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FIG. 9: (Color online) Image reconstruction for the scattering coefficient fig for an inhomogeneously 
scattering and absorbing sample and for various noise levels n. The rows show the slices x = 6h, 
13h and 20/i, where the inhomogeneities are placed. The background scattering and absorption 
coefficients are set such that flgLz = 1.6 and jla = O-ljls, the contrast in Hg varies from 2 to 3, and 
the contrast in fia varies from 2 to 5. 

using the the difference in the data functions corresponding to these measurements. On the 
other hand, the noise Sfisn in the scattering coefficient fign is estimated from Eqs. ( |T6l) and 
(in]) to be 

(\ 1/2 
i m J 

Here n labels the cell where the detected rays corresponding to the source-detection realiza- 
tion u are single scattered. Finally, the absorption coefficient is determined as the difference 
between the total attenuation coefficient and the scattering coefficient with the error 6fian 
that verifies 
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FIG. 10: (Color online) Image reconstruction for the absorption coefficient /i^ for an inhomoge- 
neously scattering and absorbing sample and for various noise levels n. The rows show the slices 
X = 6h, 13h and 20/i, where the inhomogeneities are placed. The background scattering and ab- 
sorption coefficients are set such that fisL^ = 1.6 and fia = O.l/Xs, the contrast in fig varies from 2 
to 3, and the contrast in /i^ varies from 2 to 5. 

This expresses the fact that for samples where the absorption is stronger than the scattering, 
very large noise to signal levels in the reconstructed absorption coefficient result, even for 
low noise levels in the data. This result completely explains the image reconstruction for 
the absorption coefficient presented above. The difference between Figs. [7] and [TU] is that 
the maximum value of the ratio fig/f^a is increased from 4 (for Fig. [7]) to 6 (for Fig. [TOl) . 
leading to pronounced noise in the reconstructed image. Also, although stronger scattering 
inhomogeneities are present in slice x = 13 (second row) in Fig. [TOl better image quality 
is obtained compared with the slice x = Q (ffist row), since in this case the absorbing and 
scattering inhomogeneities do not overlap and fis/f^a is smaller. The physical interpretation 
of this result is that the SSOT data function carries more signature of the scattering coeffi- 
cient than of the absorption. In other words, in SSOT, the scattering coefficient is privileged 
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FIG. 11: (Color online) Image reconstruction for the total attenuation coefficient /it for various noise 
levels n, for a data function corresponding only to single-scattered light and calculated according 
with d?]) (inverse crime). All the sample parameters are as for Fig. [HI 

as compared with the absorption coefficient. This fact originates from the RTE equation 
we employ, where the scattering coefficient has a stronger contribution than the absorption 
coefficient, and is also expressed by Eqs. ([6]) and (JTj), showing that the scattered intensity 
decays exponentially with the absorption coefficient, but has a stronger dependence on the 
scattering coefficient. In this case, besides the attenuation of the specific intensity as a result 
of absorption and scattering of photons from a given mode into other modes, there is also 
amplification of a given mode due to scattering of photons from other modes into that mode. 
We note that this dependence of the measured intensity of the scattering and absorption 
coefficients of the sample is different from the case of diffuse optical tomography. In the 
case of diffusive light propagation, the scattering of photons into and out various directional 
modes does not affect the light intensity, the diffusion equation has similar contributions 
from the scattering and absorption coefficients, and the light intensity emerging from the 
sample decays exponentially with both /i^ and /i^. 
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FIG. 12: (Color online) Image reconstruction for the scattering coefficient fig ^or various noise 
levels n, for a data function corresponding only to single-scattered light and calculated according 
with d?]) (inverse crime). All the sample parameters are as for Fig. [9l 

To verify the conclusion presented above, we perform image reconstruction for a stronger 
absorbing sample. Figs. [T4lfT6] present image reconstruction for a sample where the absorp- 
tion and scattering are spatially modulated as for Figs. IHlfTOl but the absorption coefficient 
is increased by a factor of 10, such that it becomes comparable to the scattering coeffi- 
cient. Indeed, the image reconstruction for the absorption coefficient presented in Fig. [16] 
is markedly better than in Fig. [10] and comparable to that for the scattering coefficient 
presented in Fig. [15] 

Finally, we perform image reconstruction for a stronger scattering sample, characterized 
by an optical depth of the background of figLz = 3.2 and an additional contrast in the scat- 
tering coefficient of up to 3. This is a borderline case when scattering is sufficiently strong 
so that the single-scattering approximation of SSOT may be expected to be inaccurate. The 
results for imagine reconstruction obtained for the case when the scattering and absorp- 
tion have comparable strengths are presented in Figs. [T71[T9l We obtain that even in this 
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FIG. 13: (Color online) Image reconstruction for the absorption coefficient fXa for various noise 
levels n, for a data function corresponding only to single-scattered light and calculated according 
with d?]) (inverse crime). All the sample parameters are as for Fig. 1101 

scattering regime the most relevant features in the reconstructed scattering and absorption 
coefficients remain legible. 



III. CONCLUSIONS 



We have demonstrated that the SSOT technique enables simultaneous reconstruction of 
scattering and absorption properties of mesoscopic systems. In particular, we have shown 
that while accurate, qualitative imagine reconstruction of scattering is always possible, good 
image reconstruction for absorption can be realized under the condition that scattering 
and absorption have comparable strengths. These conclusions have been reached under the 
assumption that the light propagating in the mesoscopic systems is just single scattered, but 
without making any assumption of measuring just single-scattered light. We have argued 
that better image quality for scattering as compared to absorption is possible since the SSOT 
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FIG. 14: (Color online) Image reconstruction for the total attenuation coefficient fit for an inho- 
mogeneously scattering and absorbing sample and for various noise levels n. The rows show the 
slices X = 6h, 13h and 20h, where the inhomogeneities are placed. The background scattering and 
absorption coefficients are set such that figLz = 1-6 and fia = fls, the contrast in /i^ varies from 2 
to 3, the contrast in varies from 2 to 5, and the contrast in fit varies from 2 to 4. 

data function carries a stronger signature of scattering than of absorption. 

Simultaneous reconstruction of scattering and absorption of mesoscopic systems can be 
experimentally implemented by appropriately choosing the wavelength of the illuminating 
beam, such that the effects of absorption and scattering have comparable strengths [l^ . Al- 
ternatively, the absorption characteristics of the sample can be recovered through fluorescent 
SSOT, which will be the subject of a future study. 
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FIG. 15: (Color online) Image reconstruction for the scattering coefficient fis for an inhomoge- 
neously scattering and absorbing sample and for various noise levels n. The rows show the slices 
X = 6h, 13h and 20/i, where the inhomogeneities are placed. The background scattering and ab- 
sorption coefficients are set such that figLz = 1.6 and fia = fls, the contrast in iig varies from 2 to 
3, and the contrast in //„ varies from 2 to 5. 
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FIG. 16: (Color online) Image reconstruction for the absorption coefficient /i^ for an inhomoge- 
neously scattering and absorbing sample and for various noise levels n. The rows show the slices 
X = 6h, 13h and 20/i, where the inhomogeneities are placed. The background scattering and ab- 
sorption coefficients are set such that figLz = 1.6 and fla = fis, the contrast in iig varies from 2 to 
3, and the contrast in fia varies from 2 to 5. 
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FIG. 17: (Color online) Image reconstruction for the total attenuation coefficient fit for an inho- 
mogeneously scattering and absorbing sample and for various noise levels n. The rows show the 
slices X = 6h, 13h and 20h, where the inhomogeneities are placed. The background scattering and 
absorption coefficients are set such that figLz = 3.2 and fia = fls, the contrast in /i^ varies from 2 
to 3, the contrast in varies from 2 to 5, and the contrast in fit varies from 2 to 4. 
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FIG. 18: (Color online) Image reconstruction for the scattering coefficient fis for an inhomoge- 
neously scattering and absorbing sample and for various noise levels n. The rows show the slices 
X = 6h, 13h and 20/i, where the inhomogeneities are placed. The background scattering and ab- 
sorption coefficients are set such that figLz = 3.2 and fia = fls, the contrast in iig varies from 2 to 
3, and the contrast in fia varies from 2 to 5. 
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FIG. 19: (Color online) Image reconstruction for the absorption coefficient /i^ for an inhomoge- 
neously scattering and absorbing sample and for various noise levels n. The rows show the slices 
X = 6h, 13h and 20/i, where the inhomogeneities are placed. The background scattering and ab- 
sorption coefficients are set such that figLz = 3.2 and fia = fis, the contrast in iig varies from 2 to 
3, and the contrast in fia varies from 2 to 5. 
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